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摘 要 : 太阳 高 分 辩 图 像 中 横向 速度 场 的 测量 已 AS ASAT AKA, ER 
表面 特征 的 动力 学 分 析 中 ， 但 依然 存在 测量 精度 不 够 的 问题 。 详 细 介 绍 了 
Demons 方法 ， 并 将 其 应 用 于 lm 新 真空 太阳 望远镜 (NVST) 的 高 分 辨 观测 资 
料 处 理 中 。 首 先 选 取 不 同 观 测 时 间 间 隔 和 代表 光 、 色 球 

同 波段 的 3 个 数据 集 作 为 测试 样本 ， 在 测量 到 速度 场 后 ， na 
图 像 进行 非 刚性 配 准 并 比较 结构 相似 度 (SSIM) 评价 速度 场 的 测量 精度 。 结 
KA, Demons 方法 在 小 尺度 运动 的 精细 测定 方面 ， 明 显 优 于 传统 的 FLCT 和 
并 且 采 用 光 球 和 色 球 图 像 的 亚 像 素 和 超 像 素 模拟 位 移 实验 表明 ， 

方法 的 逐 点 测量 精度 可 以 达到 0.1 像素 量 级 。 
a 太阳 高 分 辨 图 像 ; 速度 场 ; 光 流 场 ; Demons 配 准 
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SHB AIL. LAA AS HINODE P 
中 国 一 米 新 真空 太阳 望远镜 (New Vacuum Solar Telescope, {aj #K NVST) 
a 本 加 Goode Solar Telescope! Al Fit HL Swedish Solar Telescope 乌 为 代表 的 高 分 辨 
eee 7 HY Se PL 序列 资料 ， 为 太阳 物理 的 研究 带 来 重要 的 机 
遇 和 挑战 四 。 这 些 资料 直观 地 从 强度 上 展现 了 太阳 大 气 中 时 刻 发 生 的 各 种 尺度 活 
动 现象 《从 太阳 风暴 、 a 
震荡 ) 的 细节 及 变化 情况 。 高 时 空 分 辨 紊 的 图 像 不 但 能 从 形态 上 展示 太阳 大 气 
强度 变化 ， 同 时 也 为 量化 分 析 和 定量 描述 带 来 了 可 能 。 其 中 ， 通 Rae REAN 
Bn 向 速度 场 已 被 广泛 应 用 于 太阳 光 、 色 球 表面 特征 的 
B 
ARPA EDR AH EDR EY BH [Fd 1 GE A) SR (EELS BIL. OT 
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帧 之 间 亚 像素 级 或 者 几 个 像素 级 的 位 移 量 。 在 太阳 图 像 的 数据 处 理 领 域 ， 这 方面 
已 经 有 了 很 多 的 发 展 。 最 为 经 典 的 有 局 部 相关 跟踪 法 (Local Correlation 
Tracking, LCT) 四 以 及 后 来 的 基于 傅 里 叶 变 换 的 局 部 相关 跟踪 法 (Fourier Local 
Correlation Tracking, FLCT) 外 。 微 分 仿 射 速度 估计 (Differential Affine velocity 
estimator, DAVE) Ey 等 技 术 在 引入 向 场 感应 方程 后 ， 首先 被 用 于 磁 图 非 势 性 研 
究 ， 然 后 被 应 用 于 强度 像 的 测量 中 巴 币 。 此 外 文 [12-13] 提 出 了 一 种 Balltracking A 
法 ， 用 于 光 球 超 米粒 的 研究 。 虽 然 这 些 方法 都 得 到 了 一 些 很 好 的 结果 ， 但 在 测量 
o EILE NVST 这 种 像 元 分 辨 率 在 0.05arcsec, RDI K IA To 
高 分 辨 图 像 处 理 要 求 ， 因 而 有 必要 对 相关 方法 进行 更 为 深入 的 研究 。 
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移 量 得 到 。 由 于 太阳 序列 图 像 是 某 种 “ 光 " 强 度 随 着 时 间 的 变化 ， 因 此 速度 场 是 一 
种 典型 的 “ 光 流 场 "。 从 另外 一 个 方面 ， 两 帧 图 像 之 间 的 位 移 测量 也 可 以 看 作 是 一 
种 图 像 非 刚性 的 配 准 。 而 随 着 计算 机 图 像 处 理 技 术 的 发 展 ， 光 流 场 的 测量 和 图 像 
非 刚 性 配 准 技术 已 经 有 了 非常 大 的 发 展 。 
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图 像 非 刚性 配 准 方法 可 以 由 3 部 分 描述 : 〈1) 联系 源 图 像 和 目标 图 像 的 空间 
变换 ; (2) 测量 目标 图 像 和 源 图 像 相 似 性 的 相似 性 测度 G) 决定 最 优 变换 参 
数 的 优化 方法 。 

非 刚 性 配 准 方法 主要 有 基于 空间 变换 的 配 准 方法 和 基于 物理 模型 的 配 准 方法 

两 大 类 。 其 中 ， 由 文 [14] 提 出 的 基于 热力 党 中 流体 扩散 理 从 的 Demons 非 刚 性 配 
准 方法 [具有 处 理 非 光滑 变形 、 抗 噪声 性 能 好 等 优点 ， 被 广泛 应 用 于 医学 图 像 配 
HE"), HSER FFA ABA SEERA EERE 的 磁 流 体 非 均匀 局 部 变形 的 物理 特 狂 。 

本 文 用 Demons 方法 测量 NVST 高 高 分 辨 光 球 和 色 球 观测 数据 的 横 癌 速 束 度 场 。 
测量 结果 表明 该 方法 可 以 较为 精确 地 测量 太阳 高 分 辨 图 像 的 横 问 速度 场 ， 优 于 传 
统 的 FLCT 和 DAVE 方法。 

1 Demons 方法 


1.1 Demons 配 准 的 原理 和 流程 

经 典 Demons 算法 由 Thirion 提出 ， 即 “Demons-base” 算 法 中。 在 概念 上 ， 它 
和 19 世纪 Maxwell 的 实验 原理 很 相似 。 该 算法 的 基本 思想 是 假设 运动 目标 的 灰 
度 不 随时 间 改 变 ， 那 么 非 刚性 配 准 可 以 看 作 是 浮动 图 像 中 各 个 像素 问 参 考 图 像 逐 
步 扩散 的 过 程 ， 浮 动 图 像 各 个 像素 的 扩散 速度 由 参考 图 像 的 灰 度 梯度 信息 决定 。 
假设 参考 图 像 >， 浮 动 图 像 /都 是 图 像 函 数 在 某 一 时 刻 的 快照， 图 像 函数 的 灰 度 
值 保持 常数 ， 则 此 假设 可 被 模型 化 为 


在 初始 时 刻 ， 图 像 函 数 等 于 f， 经 过 一 定时 间 到 达 时 刻 后， 图 像 函 数 被 完全 变形 为 +-。 图 像 
配 准 的 过 程 就 是 要 得 出 一 个 能 驱动 了 中 的 各 个 像素 点 向 x 中 对 应 像素 点 移动 的 向 量 场 。 为 了 
得 出 驱动 力 ， 将 (1) 式 偏 微分 ， 得 到 : 

(4) 


(5) 


(4) 式 可 以 简化 为 


其 中 ， 为 从 浮动 图 像 向 参考 图 像 变 化 的 速度 场 ， 是 x 的 梯度 向 量 。 进 一 步 可 以 得 到 : 


(6) 
若 将 浮动 图 像 的 梯度 信息 作为 一 种 正 内 力 ， 参 考 图 像 的 梯度 信息 作为 负 内 力 ， 利 用 这 两 种 力 
同时 驱动 形变 ， 得 到 的 扩散 速度 为 


(7) 
ree 算法 利用 图 像 局 部 信息 来 驱动 形变 ， 在 全 局 范围 内 使 形变 连续 ， 从 而 保证 图 像 的 拓 
ane 


1.2 基于 梯度 互信 息 的 Demons 方法 
在 Demons 模型 中 ， 浮 动 图 像 各 个 像素 都 可 以 自由 移动 ， 使 得 在 浮动 图 像 中 具有 某 个 特 
灰 度 值 的 所 有 像素 有 可 能 映射 到 参考 图 像 上 的 同一 像素 点 ， 导 致 错误 的 配 准 结果 ， 为 了 使 
图 像 像素 能 沿 着 正确 的 方向 移动 , Æ (7) 式 的 基础 上 ， 采 用 梯度 互信 息 进 一 步 改进 
emons 算法 ， 在 原 扩散 速度 的 基础 上 ， 增加 了 两 幅 图 像 之 间 梯 度 互 信息 的 作用 ， 使 两 幅 
图 像 配 准 结束 的 同时 梯 ) 度 互 信息 也 达到 最 大 。 改 进 后 的 Demons 算法 扩散 速度 模型 为 
(8) 


其 中 ， 为 两 幅 图 像 的 梯度 互信 息 ; B 为 正常 数 ， 表 示 此 项 的 权重 。 
基于 图 像 灰 度 互 信息 的 配 准 方法 ， ,最 大 的 不 足 是 忽 各 了 图 像 的 空间 信息 。 就 太阳 图 像 配 
准 过 程 而 言 ， 对 相同 目标 所 成 图 像 可 能 具有 不 同 的 灰 度 ， 不 同 的 分 辨 率 ， 甚 至 图 像 本 身 大 小 
也 会 不 同 ， 但 是 对 于 同一 目标 ， 其 边界 是 确定 的 ， 不 会 发 生 很 明显 的 变化 ， 当 两 幅 图 像 配 准 
时 ， 对 应 位 置 像素 点 的 梯度 矢量 将 具有 相同 或 者 相反 的 指向 。 基 于 以 上 分 析 ， 梯 度 互信 息 采 
用 Parzen 窗 方法 进行 计算 。 


2 实验 数据 


J 


chinaXiv:201805.00538v1 


ChinaXiv 合 作 期 刊 


Im 太阳 望远镜 作为 中 国 最 大 的 太阳 望远镜 和 世界 三 大 高 分 辨 成 像 太 阳 望 远 镜 
> 一， 其 优良 的 观测 性 能 和 高 质量 的 图 像 已 经 在 很 多 方面 得 以 证 实 。 为 检验 
Demons 方法 的 可 用 性 ， 选 取 了 具有 代表 性 3 个 1m 太阳 望远镜 观测 资料 数据 集 ， 
包括 两 组 光 球 图 像 和 一 组 色 球 图 像 ， 每 组 图 像 为 前 后 两 帧 。 通 过 计算 这 两 帧 图 像 
的 光 流 场 或 者 对 它们 进行 非 刚性 配 准 ， 就 可 得 到 相应 的 速度 场 。 

数据 集 1 包括 的 两 帧 图 像 均 为 1m 太阳 望远镜 在 TiO 波段 观测 的 光 球 Level 1+ 
重建 像 ， 观 测 时 间 为 2013 年 11 月 3 日 04 时 58 分 31 秒 CUT) 和 05 时 01 分 28 
秒 (UT)〉， 时 间 间 隔 2 分 57 秒 。 视 场 大 小 为 685*650 Kt, RICA HEN 
0.04arcsec。 光 球 上 典型 的 1 km/s 的 横向 速度 将 带 来 约 6 像 元 的 位 移 。 意 味 着 整 
体 光 流 场 是 超 像素 运动 。 

数据 集 2 与 数据 集 1 是 同 目标 、 同 波段 、 同 期 的 观测 ， 只 是 第 1 帧 图 像 的 观 
测 时 间 为 05 时 1 分 10 秒 CUT) ， 与 第 2 帧 相隔 仅 18 秒 ， 大 约 是 数据 集 1 的 
1/10， 其 他 参数 与 数据 集 1 相同 。 这 样 1km/s 的 横向 速度 将 带 来 仅仅 0.6 像 元 的 位 
移 ， 即 光 流 场 为 亚 像素 运动 。 

数据 集 3 的 两 帧 图 像 为 Im 太阳 望远镜 在 了 .波段 的 Level 1+ 重 建 像 。 观 测 时 
间 为 2017 年 8 月 6 日 09 时 59 分 $3 秒 和 10 时 00 分 25 秒 ， 时 间 间 隔 为 32 秒 。 视 
场 大 小 为 238*225 像 元 ， 像 元 分 辨 率 0.13arcsec。 色 球 上 5km/s 的 横向 速度 将 带 


图 1 左 为 数据 集 1 和 2 中 的 第 二 帧 ， 右 为 数据 集 3 中 的 第 二 帧 
Fig. 1 The second frame of dataset 1 and 2 (left) and the second frame of dataset 3 (right) 


从 图 1 可 以 看 到 ， 光 球 像 数 据 集 1 和 2 是 一 个 太阳 活动 区 的 局 部 ， 特 征 结构 履 


盖 了 米粒 、 黑 子 本 影 和 半 影 。 数 据 集 3 是 典型 的 色 球 纤 维 ， 图 像 结 构 本 吴 就 有 很 
强 的 自 相似 性 。 
3 ”数据 处 理 结果 


对 上 述 3 个 数据 集 采 用 Demons 方法 计算 第 2 帧 相对 于 第 1 帧 的 逐 像素 位 移 量 ， 
即 光 流 场 ， 然 后 根据 相应 的 时 间 间 隔 归 算 到 逐 点 的 横向 速度 场 。 图 2 显示 了 数据 
集 1 中 一 个 小 区 域 的 两 帧 图 以 及 用 Demons 方法 计算 出 的 位 移 矢 量 。 米 粒 的 细微 
运动 可 以 清晰 地 反映 在 图 上 。 


图 2 数据 集 1 中 同一 区 域 的 变化 及 用 Demons 方法 测量 到 的 矢量 场 
Fig. 2 The same region of two frames in dataset | and the vector fields measured by Demons 
为 了 测试 Demons 的 结果 并 与 目前 太阳 观测 资料 处 理 中 常用 的 FLCT 和 
DAVE 方法 进行 比 对 ， 在 分 别 采 用 不 同方 法 得 到 光 流 场 后 ， 采 用 双 线 性 插值 将 第 


2 帧 与 第 1 帧 进行 非 刚 性 配 准 ， 并 利用 图 像 结构 相似 度 指 数 〈Structural Similarity 
Index Measurement, SSIM) 09 测试 配 准 的 精度 。 较 高 的 配 准 精度 即 意味 着 较 高 的 
光 流 场 〈 速 度 场 ) 测量 精度 。 

对 于 x，y 两 张 图 像 ， 其 图 像 结构 相似 度 指 数 由 下 列 公 式 求 得 : 


其 中 ， 和 为 两 张 图 像 的 均值 ， 和 为 各 自 的 方差 ， 为 彼此 的 协 方 差 。Cl 和 CG; 是 
为 了 防止 分 母 为 0 而 设置 的 系数 ， 对 于 强度 归 一 化 的 图 像 其 值 分 别 为 0.01 和 
0.03。 图 像 结构 相似 度 指 数 的 取 值 范围 为 0 到 1， 相 似 程度 越 大 ， 其 值 越 大 。 如 
果 两 张 图 像 完 全 相同 ， 则 为 1。 

为 了 逐 像素 比 对 非 刚 性 配 准 后 的 两 帧 图 像 的 相似 度 ， 采 用 11*11 像素 邻 域 的 
窗口 平均 作为 该 点 的 局 部 图 像 结构 相似 度 指数 值 。 对 于 配 准 的 结果 ， 如 果 有 更 多 
点 的 图 像 结构 相似 度 指数 值 接近 1， 说 明 该 方法 总 体位 移 测量 更 为 准确 ， 配 准 更 
好 。 而 图 像 中 所 有 像素 点 图 像 结构 相似 度 指数 值 的 平均 作为 两 帧 图 像 的 平均 相似 
度 值 ， 也 代表 整体 测量 水 平 。 

表 1 给 出 了 3 个 数据 集 在 未 配 准 前 和 使 用 FLCT、DAVE 和 Demons 方法 分 别 
测量 并 非 刚 性 配 准 后 与 第 一 帧 图 像 比较 的 平均 图 像 结构 相似 度 指数 值 。 


表 1 不 同方 法 配 准 后 的 平均 SSIM 值 
Tab. 1 The SSIM index after non-rigid registration with different method 
数据 集 1 数据 集 2 数据 集 3 


未 配 准 0. 7410 0. 9884 0. 6731 
FLCT 0. 7440 0. 9885 0. 6881 
DAVE 0. 8332 0.9921 0. 7532 
Demons 0.9072 0. 9965 0. 8298 


从 表 1 可 以 看 出 ， 原 始 数据 彼此 的 相似 度 不 但 取决 于 时 间 间 陋 ， 也 取决 于 特 
征 结构 的 运动 速度 。 对 于 色 球 数据 ， 快 速 的 运动 导致 在 半分 钟 内 图 像 就 有 较 大 的 
变化 。 另 外 ， 使 用 FLCT 配 准 后 ， 图 像 结构 相似 度 指数 提升 非常 有 限 ， 即 只 测量 
ee new DAVE 方法 的 测量 结果 有 显著 的 提升 ， 但 Demons 方法 更 
ve 


数据 集 1 间隔 时 间 较 长 ， 而 数据 集 3 的 色 球 运动 较 快 。 在 这 两 个 数据 集中 ， 
特征 结构 除了 局 部 的 运动 外 ， 还 有 新 结构 的 浮现 和 旧 结 构 的 消失 ， 因 此 图 像 结 构 
相似 度 指 数 相 对 偏 小 。 而 数据 集 2 中 两 帧 图 像 时 间 间 隔 仅 为 18 s， 且 光 球 运动 较 
为 绥 慢 ， 图 像 差 异 微 小 ， 图 像 结构 相似 度 指数 也 明显 高 于 其 他 两 个 数据 集 。 但 无 
论 原 始 图 像 情 况 如 何 ，Demons 方法 对 于 不 同 的 数据 集 均 表 现 出 良好 的 处 理 能 力 。 

图 3 显示 了 对 数据 集 1 应 用 3 种 方法 配 准 后 的 逐 点 图 像 结构 相似 度 指数 的 直 
方 图 概率 分 布 。 可 以 看 到 Demons 方法 配 准 后 图 像 结 构 相 似 度 指数 值 的 整体 分 布 
更 为 靠近 1， 明 显 好 于 DAVE 和 FLCT。 同 时 FLCT 的 结果 和 未 配 准 图 像 是 非常 
EN: 即 测量 效果 并 不 显著 。 这 与 表 1 的 结果 是 一 致 的 。 数 据 集 2 和 3 的 图 像 

结构 相似 度 指数 概率 分 布 也 反映 出 同样 的 趋势 。 
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图 3 数据 集 1 的 不 同方 法 配 准 后 的 逐 点 SSIM 概率 分 布 
Fig. 3 The probability distribution of SSIM index with different registration method applying dataset 1 


从 最 终 速度 场 的 测量 情况 看 ，3 个 数据 集 都 呈现 出 同样 的 趋势 : FLCT 基本 
上 只 能 在 位 移 较 大 的 情况 下 有 测量 值 ， 而 DAVE 的 测量 值 系统 性 小 于 Demons 的 
测量 值 ， 大 多 数 情况 只 有 Demons 方法 测量 精度 的 一 半 。 


4 ”模拟 数据 测试 


Demons 方法 虽然 结果 好 于 DAVE 和 FLCT， 但 第 2 帧 经 过 非 刚 性 配 准 后 也 并 
没有 做 到 完全 的 匹配 ， 即 图 像 结构 相似 度 指 数 并 不 等 于 1。 其 中 最 为 重要 的 原因 
无 论 是 光 球 还 是 色 球 ， 经 过 一 段 时 间 间 隔 后 ， 不 但 原 有 的 米粒 、 黑 子 或 者 色 球 纤 
ctrl 同时 也 有 新 特征 的 产生 和 旧 结 构 的 消失 ， 这 就 不 是 仅 从 光 流 场 能 

计算 的 。 


为 了 进一步 比较 Demons 和 DAVE 的 光 流 场 测量 精度 ， 测 试 他 们 在 不 同 太 阳 
RR ERMER) 上 亚 像 素 和 超 像 素 位 移 的 可 靠 性 ， 设 计 了 两 个 模拟 实验 ， 即 
分 别 将 数据 集 1 和 数据 集 3 中 的 第 1 帧 图 像 在 式 和 了 方向 整体 移动 一 个 固定 量 ， 
从 而 产生 第 2 帧 刚性 运动 的 图 像 ， 然 后 分 别 再 用 两 个 方法 进行 测量 其 与 第 一 帧 图 
像 的 光 流 场 。 由 于 模拟 的 是 全 局 性 位 移 ， 因 此 其 逐 点 光 流 场 的 均值 和 标准 差 即 可 
以 很 好 地 反映 测量 的 精度 。 

R 2 给 出 了 整体 位 移 为 (5.7, -0.2) 时 ， 这 两 个 方法 在 对 和 了 方向 上 测量 的 统 
计 结 


表 2 模拟 数据 使 用 DAVE 和 Demons 方法 测量 的 统计 结果 
Tab. 2 The statistics results of DAVE and Demons for simulation data 


HIRA A 光 球 了 方向 色 球 X 方 向 色 球 了 方向 
DAVE 均值 5. 1005 -0. 4316 1. 7071 -0. 4550 
DAVE 标准 差 1. 3301 0. 9575 1. 1969 1. 9099 
Demons 均值 5. 6946 -0. 2075 5. 6945 -0. 2069 
Demons 标准 差 0. 0856 0. 1048 0. 0716 0. 1694 


可 以 看 到 ， 与 预 设 值 相 比 ，Demons 的 均值 和 标准 差 表 现 都 明显 好 于 DAVE 
方法 ， 逐 点 测量 精度 在 0.1 像素 。 而 DAVE 不 但 标准 差 比 Demons 方法 高 一 个 量 
级 ， 而 且 均 值 也 有 较 大 的 系统 偏差 ， 尤 其 对 于 自 相 关 较 强 的 色 球 图 像 。 一 般 认为 
DAVE 在 做 超 像 素 测量 时 误差 较 大 ， 但 测试 结果 表明 DAVE 受到 图 像 本 身 的 特征 
结构 的 影响 也 是 很 大 的 。 这 与 DAVE 本 质 上 也 是 使 用 相关 方法 来 测量 有 关 。 
FLCT 由 于 能 计算 的 像素 太 少 ， 缺 乏 统计 意义 ， 因 此 被 排除 在 外 。 

同时 ， 由 于 我 们 模拟 整 图 的 刚性 位 移 ， 因 此 也 可 以 用 标准 的 整 幅 图 像 互 相关 
方法 直接 测量 两 组 图 像 的 刚性 位 移 。 光 球 像 的 测量 结果 为 (5.7578, -0.1814) ， 


而 色 球 像 的 结果 为 (5.7446，-0.0654)， 也 比 用 Demons 方法 逐 点 测量 后 的 均值 结 
果 明 显 偏差 要 大 。 这 意味 着 ，Demons 方法 不 但 适用 于 非 刚 性 的 配 准 ， 同 时 也 可 
以 利用 均值 测量 图 像 刚 性 的 位 移 ， 用 于 图 像 的 刚性 配 准 。 尤 其 对 于 自 相 似 性 较 强 
的 目标 ， 如 色 球 纤维 、 太 阳 边 缘 等 ， 这 一 方法 有 较 强 的 优势 。 


5 结论 


本 文 详细 描述 了 Demons 方法 ， 并 利用 其 测量 不 同 目标 和 不 同位 移 量 的 太阳 
高 分 辨 像 的 横 疝 速度 场 。 对 测量 结果 进行 图 像 非 刚性 配 准 后 ， 图 像 结 构 相 似 度 指 
数 有 了 明显 的 提高 ， 从 而 说 明 位 移 的 测量 是 更 为 准确 ， 表 明了 这 一 方法 的 可 行 性 。 
通过 和 DAVE 和 FLCT 方法 进行 对 比 ，Demons 方法 表现 出 较为 明显 的 优势 。 而 
光 球 和 色 球 数据 的 亚 像 素 和 超 像 素 模拟 位 移 实验 表明 ，Demons 逐 点 测量 精度 在 
0.1 像素 量 级 。 这 意味 着 在 Im 太阳 望远镜 18s 间隔 的 光 球 像 上 可 以 得 到 200 m/s 
青 度 的 横向 速度 ， 而 在 32s 的 色 球 像 上 可 以 得 到 500 m/s 的 速度 精度 。 另 外 ， 对 
于 图 像 结构 自 相关 性 较 大 的 图 像 配 准 ， 这 一 方法 无 论 在 非 刚性 和 刚性 配 准 方面 相 
对 于 传统 互相 关 方 法 有 明显 的 有 优势 。 

对 比 研究 同时 也 表明 ，FLCT 基本 不 适合 测量 高 分 辨 像 ， 而 DAVE 直接 用 于 
强度 像 的 测量 也 难以 得 到 精确 的 结果 。 

Demons 方法 在 亚 像素 和 超 像素 测量 上 都 能 获得 较 高 的 精度 ， 但 是 由 于 其 算 
法 的 复杂 性 ， 使 得 计算 机 开销 较 大 ， 计 算 时 间 较 长 ， 并 不 太 适 合 大 尺寸 高 分 辨 图 
像 的 整体 运算 。 同 时 处 理 的 图 像 有 较为 明显 的 边界 效应 ， 因 此 在 实际 处 理 的 时 候 ， 
需要 对 关注 的 区 域 做 一 定 的 边缘 扩展 。 

必须 指出 的 是 ， 这 种 利用 光 流 场 来 测量 横向 速度 场 从 本 质 上 是 假设 物质 的 运 
动 表现 为 图 像 强 度 的 横向 变化 ， 而 且 这 也 仅仅 是 在 观测 平面 上 的 投影 ， 并 不 是 真 
正 的 物质 流动 方向 。 同 时 如 果 有 新 图 像 结 构 的 产生 和 旧 结 构 的 消失 ， 从 本 质 上 这 
种 测量 方式 是 有 缺陷 的 ， 因 此 需要 高 时 间 分 辨 率 的 观测 ， 以 减少 这 种 变化 。 但 更 
短 的 时 间 间 隔 意味 着 更 小 的 位 移 量 。 

理论 上 ， 这 一 方法 可 以 用 于 处 理 其 他 望远镜 的 高 分 辨 观测 图 像 ， 也 可 以 用 于 
一 些 低 分 辨 数据 ， 或 者 其 他 类 型 的 图 像 上 ， 这 些 还 需要 进一步 的 尝试 。 
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Measurement of transverse velocity field of NVST 


solar high resolution image 
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Abstract: The measurement of the transverse velocity field in the high resolution 
image of the Sun has been widely used in the dynamic analysis of the surface 
characteristics of the Chromosphere and Photosphere, however, the problem of 
insufficient measurement accuracy still exists. 

In this paper, the Demons method is introduced in detail and applied to the high 
resolution observation data processing of one meter vacuum solar telescope (NVST) 
in China. First, we select 
three data sets representing different bands of light and color balls at different 
observational intervals as test samples. When the velocity field is measured, the 
measurement accuracy of the velocity field is evaluated by the non-rigid registration 
with the image of previous time and the comparison of the structural similarity 
(SSIM). Experimental results show that Demons method is superior to traditional 
FLCT and Dave methods in the fine measurement of small scale motion, and the 
experimental results of sub-pixel and hyper-pixel analog displacements using 
Photosphere and Chromosphere image show that the point measurement accuracy can 
reach 0.1 pixel magnitude by our method. 

Key words: Solar high-resolution image; Velocity field; Optical flow field ; Demons 
registration 


